The AccelerAge framework: a new statistical approach to predict biological age based on time-to-event data.

Aging is a multifaceted and intricate physiological process characterized by a gradual decline in functional capacity, leading to increased susceptibility to diseases and mortality. While chronological age serves as a strong risk factor for age-related health conditions, considerable heterogeneity exists in the aging trajectories of individuals, suggesting that biological age may provide a more nuanced understanding of the aging process. However, the concept of biological age lacks a clear operationalization, leading to the development of various biological age predictors without a solid statistical foundation. This paper addresses these limitations by proposing a comprehensive operationalization of biological age, introducing the "AccelerAge" framework for predicting biological age, and introducing previously underutilized evaluation measures for assessing the performance of biological age predictors. The AccelerAge framework, based on Accelerated Failure Time (AFT) models, directly models the effect of candidate predictors of aging on an individual's survival time, aligning with the prevalent metaphor of aging as a clock. We compare predictors based on the AccelerAge framework to a predictor based on the GrimAge predictor, which is considered one of the best-performing biological age predictors, using simulated data as well as data from the UK Biobank and the Leiden Longevity Study. Our approach seeks to establish a robust statistical foundation for biological age clocks, enabling a more accurate and interpretable assessment of an individual's aging status.

Epigenome-wide association study of dietary fatty acid intake.

BACKGROUND: Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS). RESULTS: DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively). CONCLUSIONS: Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

Identified senescence endotypes in aged cartilage are reflected in the blood metabolome.

Heterogeneous accumulation of senescent cells expressing the senescence-associated secretory phenotype (SASP) affects tissue homeostasis which leads to diseases, such as osteoarthritis (OA). In this study, we set out to characterize heterogeneity of cellular senescence within aged articular cartilage and explored the presence of corresponding metabolic profiles in blood that could function as representative biomarkers. Hereto, we set out to perform cluster analyses, using a gene-set of 131 senescence genes (N = 57) in a previously established RNA sequencing dataset of aged articular cartilage and a generated metabolic dataset in overlapping blood samples. Using unsupervised hierarchical clustering and pathway analysis, we identified two robust cellular senescent endotypes. Endotype-1 was enriched for cell proliferating pathways, expressing forkhead box protein O4 (FOXO4), RB transcriptional corepressor like 2 (RBL2), and cyclin-dependent kinase inhibitor 1B (CDKN1B); the FOXO mediated cell cycle was identified as possible target for endotype-1 patients. Endotype-2 showed enriched inflammation-associated pathways, expressed by interleukin 6 (IL6), matrix metallopeptidase (MMP)1/3, and vascular endothelial growth factor (VEGF)C and SASP pathways were identified as possible targets for endotype-2 patients. Notably, plasma-based metabolic profiles in overlapping blood samples (N = 21) showed two corresponding metabolic clusters in blood. These non-invasive metabolic profiles could function as biomarkers for patient-tailored targeting of senescence in OA.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

Circulating cell-free DNA in health and disease - the relationship to health behaviours, ageing phenotypes and metabolomics.

Circulating cell-free DNA (cf-DNA) has emerged as a promising biomarker of ageing, tissue damage and cellular stress. However, less is known about health behaviours, ageing phenotypes and metabolic processes that lead to elevated cf-DNA levels. We sought to analyse the relationship of circulating cf-DNA level to age, sex, smoking, physical activity, vegetable consumption, ageing phenotypes (physical functioning, the number of diseases, frailty) and an extensive panel of biomarkers including blood and urine metabolites and inflammatory markers in three human cohorts (N = 5385; 17-82 years). The relationships were assessed using correlation statistics, and linear and penalised regressions (the Lasso), also stratified by sex.cf-DNA levels were significantly higher in men than in women, and especially in middle-aged men and women who smoke, and in older more frail individuals. Correlation statistics of biomarker data showed that cf-DNA level was higher with elevated inflammation (C-reactive protein, interleukin-6), and higher levels of homocysteine, and proportion of red blood cells and lower levels of ascorbic acid. Inflammation (C-reactive protein, glycoprotein acetylation), amino acids (isoleucine, leucine, tyrosine), and ketogenesis (3-hydroxybutyrate) were included in the cf-DNA level-related biomarker profiles in at least two of the cohorts.In conclusion, circulating cf-DNA level is different by sex, and related to health behaviour, health decline and metabolic processes common in health and disease. These results can inform future studies where epidemiological and biological pathways of cf-DNA are to be analysed in details, and for studies evaluating cf-DNA as a potential clinical marker.

Metabolomic profiles predict individual multidisease outcomes.

Risk stratification is critical for the early identification of high-risk individuals and disease prevention. Here we explored the potential of nuclear magnetic resonance (NMR) spectroscopy-derived metabolomic profiles to inform on multidisease risk beyond conventional clinical predictors for the onset of 24 common conditions, including metabolic, vascular, respiratory, musculoskeletal and neurological diseases and cancers. Specifically, we trained a neural network to learn disease-specific metabolomic states from 168 circulating metabolic markers measured in 117,981 participants with ~1.4 million person-years of follow-up from the UK Biobank and validated the model in four independent cohorts. We found metabolomic states to be associated with incident event rates in all the investigated conditions, except breast cancer. For 10-year outcome prediction for 15 endpoints, with and without established metabolic contribution, a combination of age and sex and the metabolomic state equaled or outperformed established predictors. Moreover, metabolomic state added predictive information over comprehensive clinical variables for eight common diseases, including type 2 diabetes, dementia and heart failure. Decision curve analyses showed that predictive improvements translated into clinical utility for a wide range of potential decision thresholds. Taken together, our study demonstrates both the potential and limitations of NMR-derived metabolomic profiles as a multidisease assay to inform on the risk of many common diseases simultaneously.

Age, Sex, and BMI Influence on Copper, Zinc, and Their Major Serum Carrier Proteins in a Large European Population Including Nonagenarian Offspring From MARK-AGE Study.

The analysis of copper (Cu) and zinc (Zn) along with their major serum carriers, albumin (Alb) and ceruloplasmin (Cp), could provide information on the capacity of humans to maintain homeostasis of metals (metallostasis). However, their relationship with aging, sex, body mass index, as well as with nutritional and inflammatory markers was never investigated in a large-scale study. Here, we report results from the European large-scale cross-sectional study MARK-AGE in which Cu, Zn, Alb, Cp, as well as nutritional and inflammatory parameters were determined in 2424 age-stratified participants (35-75 years), including the general population (RASIG), nonagenarian offspring (GO), a well-studied genetic model of longevity, and spouses of GO (SGO). In RASIG, Cu to Zn ratio and Cp to Alb ratio were higher in women than in men. Both ratios increased with aging because Cu and Cp increased and Alb and Zn decreased. Cu, Zn, Alb, and Cp were found associated with several inflammatory as well as nutritional biomarkers. GO showed higher Zn levels and higher Zn to Alb ratio compared to RASIG, but we did not observe significant differences with SGO, likely as a consequence of the low sample size of SGO and the shared environment. Our results show that aging, sex, body mass index, and GO status are characterized by different levels of Cu, Zn, and their serum carrier proteins. These data and their relationship with inflammatory biomarkers support the concept that loss of metallostasis is a characteristic of inflammaging.

Improved selection of participants in genetic longevity studies: family scores revisited.

BACKGROUND: Although human longevity tends to cluster within families, genetic studies on longevity have had limited success in identifying longevity loci. One of the main causes of this limited success is the selection of participants. Studies generally include sporadically long-lived individuals, i.e. individuals with the longevity phenotype but without a genetic predisposition for longevity. The inclusion of these individuals causes phenotype heterogeneity which results in power reduction and bias. A way to avoid sporadically long-lived individuals and reduce sample heterogeneity is to include family history of longevity as selection criterion using a longevity family score. A main challenge when developing family scores are the large differences in family size, because of real differences in sibship sizes or because of missing data. METHODS: We discussed the statistical properties of two existing longevity family scores: the Family Longevity Selection Score (FLoSS) and the Longevity Relatives Count (LRC) score and we evaluated their performance dealing with differential family size. We proposed a new longevity family score, the mLRC score, an extension of the LRC based on random effects modeling, which is robust for family size and missing values. The performance of the new mLRC as selection tool was evaluated in an intensive simulation study and illustrated in a large real dataset, the Historical Sample of the Netherlands (HSN). RESULTS: Empirical scores such as the FLOSS and LRC cannot properly deal with differential family size and missing data. Our simulation study showed that mLRC is not affected by family size and provides more accurate selections of long-lived families. The analysis of 1105 sibships of the Historical Sample of the Netherlands showed that the selection of long-lived individuals based on the mLRC score predicts excess survival in the validation set better than the selection based on the LRC score . CONCLUSIONS: Model-based score systems such as the mLRC score help to reduce heterogeneity in the selection of long-lived families. The power of future studies into the genetics of longevity can likely be improved and their bias reduced, by selecting long-lived cases using the mLRC.

Broad phenotype of cysteine-altering NOTCH3 variants in UK Biobank: CADASIL to nonpenetrance.

OBJECTIVE: To determine the small vessel disease spectrum associated with cysteine-altering NOTCH3 variants in community-dwelling individuals by analyzing the clinical and neuroimaging features of UK Biobank participants harboring such variants. METHODS: The exome and genome sequencing datasets of the UK Biobank (n = 50,000) and cohorts of cognitively healthy elderly (n = 751) were queried for cysteine-altering NOTCH3 variants. Brain MRIs of individuals harboring such variants were scored according to Standards for Reporting Vascular Changes on Neuroimaging criteria, and clinical information was extracted with ICD-10 codes. Clinical and neuroimaging data were compared to age- and sex-matched UK Biobank controls and clinically diagnosed patients from the Dutch cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) registry. RESULTS: We identified 108 individuals harboring a cysteine-altering NOTCH3 variant (2.2 of 1,000), of whom 75% have a variant that has previously been reported in CADASIL pedigrees. Almost all variants were located in 1 of the NOTCH3 protein epidermal growth factor-like repeat domains 7 to 34. White matter hyperintensity lesion load was higher in individuals with NOTCH3 variants than in controls (p = 0.006) but lower than in patients with CADASIL with the same variants (p < 0.001). Almost half of the 24 individuals with brain MRI had a Fazekas score of 0 or 1 up to age 70 years. There was no increased risk of stroke. CONCLUSIONS: Although community-dwelling individuals harboring a cysteine-altering NOTCH3 variant have a higher small vessel disease MRI burden than controls, almost half have no MRI abnormalities up to age 70 years. This shows that NOTCH3 cysteine altering variants are associated with an extremely broad phenotypic spectrum, ranging from CADASIL to nonpenetrance.

Systematic Evaluation of Normalization Methods for Glycomics Data Based on Performance of Network Inference.

Glycomics measurements, like all other high-throughput technologies, are subject to technical variation due to fluctuations in the experimental conditions. The removal of this non-biological signal from the data is referred to as normalization. Contrary to other omics data types, a systematic evaluation of normalization options for glycomics data has not been published so far. In this paper, we assess the quality of different normalization strategies for glycomics data with an innovative approach. It has been shown previously that Gaussian Graphical Models (GGMs) inferred from glycomics data are able to identify enzymatic steps in the glycan synthesis pathways in a data-driven fashion. Based on this finding, here, we quantify the quality of a given normalization method according to how well a GGM inferred from the respective normalized data reconstructs known synthesis reactions in the glycosylation pathway. The method therefore exploits a biological measure of goodness. We analyzed 23 different normalization combinations applied to six large-scale glycomics cohorts across three experimental platforms: Liquid Chromatography - ElectroSpray Ionization - Mass Spectrometry (LC-ESI-MS), Ultra High Performance Liquid Chromatography with Fluorescence Detection (UHPLC-FLD), and Matrix Assisted Laser Desorption Ionization - Furier Transform Ion Cyclotron Resonance - Mass Spectrometry (MALDI-FTICR-MS). Based on our results, we recommend normalizing glycan data using the 'Probabilistic Quotient' method followed by log-transformation, irrespective of the measurement platform. This recommendation is further supported by an additional analysis, where we ranked normalization methods based on their statistical associations with age, a factor known to associate with glycomics measurements.

Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases.

Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (N = 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of IKZF1 decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.

Prevalence and Loads of Torquetenovirus in the European MARK-AGE Study Population.

Torquetenovirus (TTV) viremia has been associated with increased mortality risk in the elderly population. This work aims to investigate TTV viremia as a potential biomarker of immunosenescence. We compared levels of circulating TTV in 1813 participants of the MARK-AGE project, including human models of delayed (offspring of centenarians [GO]) and premature (Down syndrome [DS]) immunosenescence. The TTV load was positively associated with age, cytomegalovirus (CMV) antibody levels, and the Cu/Zn ratio and negatively associated with platelets, total cholesterol, and total IgM. TTV viremia was highest in DS and lowest in GO, with intermediate levels in the SGO (spouses of GO) and RASIG (Randomly Recruited Age-Stratified Individuals From The General Population) populations. In the RASIG population, TTV DNA loads showed a slight negative association with CD3+T-cells and CD4+T-cells. Finally, males with ≥4log TTV copies/mL had a higher risk of having a CD4/CD8 ratio<1 than those with lower viremia (odds ratio [OR] = 2.85, 95% confidence interval [CI]: 1.06-7.62), as well as reduced CD3+ and CD4+T-cells compared to males with lower replication rates (<4log), even after adjusting for CMV infection. In summary, differences in immune system preservation are reflected in the models of delayed and premature immunosenescence, displaying the best and worst control over TTV replication, respectively. In the general population, TTV loads were negatively associated with CD4+ cell counts, with an increased predisposition for an inverted CD4/CD8 ratio for individuals with TTV loads ≥4log copies/mL, thus promoting an immune risk phenotype.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals.

BACKGROUND: Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans. OBJECTIVE: The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes. METHODS: A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes. RESULTS: The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs. CONCLUSIONS: We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.

Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage.

OBJECTIVE: Multiple single-nucleotide polymorphisms (SNPs) conferring susceptibility to osteoarthritis (OA) mark imbalanced expression of positional genes in articular cartilage, reflected by unequally expressed alleles among heterozygotes (allelic imbalance [AI]). We undertook this study to explore the articular cartilage transcriptome from OA patients for AI events to identify putative disease-driving genetic variation. METHODS: AI was assessed in 42 preserved and 5 lesioned OA cartilage samples (from the Research Arthritis and Articular Cartilage study) for which RNA sequencing data were available. The count fraction of the alternative alleles among the alternative and reference alleles together (φ) was determined for heterozygous individuals. A meta-analysis was performed to generate a meta-φ and P value for each SNP with a false discovery rate (FDR) correction for multiple comparisons. To further validate AI events, we explored them as a function of multiple additional OA features. RESULTS: We observed a total of 2,070 SNPs that consistently marked AI of 1,031 unique genes in articular cartilage. Of these genes, 49 were found to be significantly differentially expressed (fold change <0.5 or >2, FDR <0.05) between preserved and paired lesioned cartilage, and 18 had previously been reported to confer susceptibility to OA and/or related phenotypes. Moreover, we identified notable highly significant AI SNPs in the CRLF1, WWP2, and RPS3 genes that were related to multiple OA features. CONCLUSION: We present a framework and resulting data set for researchers in the OA research field to probe for disease-relevant genetic variation that affects gene expression in pivotal disease-affected tissue. This likely includes putative novel compelling OA risk genes such as CRLF1, WWP2, and RPS3.

Sex Differences in Genetic Associations With Longevity.

IMPORTANCE: Sex differences in genetic associations with human longevity remain largely unknown; investigations on this topic are important for individualized health care. OBJECTIVE: To explore sex differences in genetic associations with longevity. DESIGN SETTING AND PARTICIPANTS: This population-based case-control study used sex-specific genome-wide association study and polygenic risk score (PRS) analyses to examine sex differences in genetic associations with longevity. Five hundred sixty-four male and 1614 female participants older than 100 years were compared with a control group of 773 male and 1526 female individuals aged 40 to 64 years. All were Chinese Longitudinal Healthy Longevity Study participants with Han ethnicity who were recruited in 1998 and 2008 to 2014. MAIN OUTCOMES AND MEASURES: Sex-specific loci and pathways associated with longevity and PRS measures of joint effects of sex-specific loci. RESULTS: Eleven male-specific and 11 female-specific longevity loci (P < 10-5) and 35 male-specific and 25 female-specific longevity loci (10-5 ≤ P < 10-4) were identified. Each of these loci's associations with longevity were replicated in north and south regions of China in one sex but were not significant in the other sex (P = .13-.97), and loci-sex interaction effects were significant (P < .05). The associations of loci rs60210535 of the LINC00871 gene with longevity were replicated in Chinese women (P = 9.0 × 10-5) and US women (P = 4.6 × 10-5) but not significant in Chinese and US men. The associations of the loci rs2622624 of the ABCG2 gene were replicated in Chinese women (P = 6.8 × 10-5) and European women (P = .003) but not significant in both Chinese and European men. Eleven male-specific pathways (inflammation and immunity genes) and 34 female-specific pathways (tryptophan metabolism and PGC-1α induced) were significantly associated with longevity (P < .005; false discovery rate < 0.05). The PRS analyses demonstrated that sex-specific associations with longevity of the 4 exclusive groups of 11 male-specific and 11 female-specific loci (P < 10-5) and 35 male-specific and 25 female-specific loci (10-5 ≤P < 10-4) were jointly replicated across north and south discovery and target samples. Analyses using the combined data set of north and south showed that these 4 groups of sex-specific loci were jointly and significantly associated with longevity in one sex (P = 2.9 × 10-70 to 1.3 × 10-39) but not jointly significant in the other sex (P = .11 to .70), while interaction effects between PRS and sex were significant (P = 4.8 × 10-50 to 1.2 × 10-16). CONCLUSION AND RELEVANCE: The sex differences in genetic associations with longevity are remarkable, but have been overlooked by previously published genome-wide association studies on longevity. This study contributes to filling this research gap and provides a scientific basis for further investigating effects of sex-specific genetic variants and their interactions with environment on healthy aging, which may substantially contribute to more effective and targeted individualized health care for male and female elderly individuals.

Survival analysis with delayed entry in selected families with application to human longevity.

In the field of aging research, family-based sampling study designs are commonly used to study the lifespans of long-lived family members. However, the specific sampling procedure should be carefully taken into account in order to avoid biases. This work is motivated by the Leiden Longevity Study, a family-based cohort of long-lived siblings. Families were invited to participate in the study if at least two siblings were 'long-lived', where 'long-lived' meant being older than 89 years for men or older than 91 years for women. As a result, more than 400 families were included in the study and followed for around 10 years. For estimation of marker-specific survival probabilities and correlations among life times of family members, delayed entry due to outcome-dependent sampling mechanisms has to be taken into account. We consider shared frailty models to model left-truncated correlated survival data. The treatment of left truncation in shared frailty models is still an open issue and the literature on this topic is scarce. We show that the current approaches provide, in general, biased estimates and we propose a new method to tackle this selection problem by applying a correction on the likelihood estimation by means of inverse probability weighting at the family level.

Zinc-Induced Metallothionein in Centenarian Offspring From a Large European Population: The MARK-AGE Project.

Metallothionein (MT) family are cysteine-rich proteins that regulate zinc (Zn) homeostasis and protect against oxidative damage. Studies in transgenic mice have shown that MT favorably influence longevity, although their role in human aging is not completely understood. Within the European multicenter study MARK-AGE, we analyzed MT induction after Zn treatment in peripheral blood mononuclear cells (PBMCs) and its relation with redox biomarkers in 2,936 age-stratified subjects (35-75 years) including the general population (RASIG), centenarian offspring (GO), and their spouses (SGO). We found that the lymphocyte capability to induce MT in response to Zn is not affected by aging. However, GO participants showed lower Zn-induced MT and increased basal expression of MT1A, MT1X, and ZnT-1 genes than RASIG subjects. Moreover, Zn-induced MT levels were found to be inversely related with oxidative stress markers (plasma protein carbonyls, 3-nitrotyrosine, and malondialdehyde) in the whole population, but not in GO subjects. In conclusion, our results support the hypothesis that the response to Zn is attenuated in PBMCs of centenarian offspring compared to the general population as a consequence of a tighter control of Zn homeostasis which is likely to provide them constant protection against stress stimuli over the whole lifespan.

Automated quantification of metabolites in blood-derived samples by NMR.

NMR is widely applied in the field of metabolomics due to the quantitative nature of the technology and the reproducible data generated. However, because of severe spectral crowding, quantifying individual metabolites in body fluids such as serum and plasma remains a challenge. In this study, a method to automatically annotate and quantify a number of small metabolites in human serum and EDTA plasma is introduced. It combines the superior signal-to-noise ratio of the commonly applied CPMG and NOESY1D pulse sequences with the superior resolution of the 2D JRES experiment to construct a model that extracts the metabolite concentrations directly from the 1D spectra without tedious deconvolution. The performance of the method was assessed by comparing the calculated areas of the various glucose peaks with known clinical values, by comparing several peaks of the same metabolite (extracted versus non-extracted), and by comparing areas obtained from various NMR pulse sequences. Additionally, the models were tested on independent datasets. It was found that for many metabolites peaks could be assigned that show a consistent behavior, indicating a precise quantification. The same method should be applicable to other biofluids with a stable composition and pH, such as CSF fluid, cell extracts, and cell media.

Secondary phenotype analysis in ascertained family designs: application to the Leiden longevity study.

The case-control design is often used to test associations between the case-control status and genetic variants. In addition to this primary phenotype, a number of additional traits, known as secondary phenotypes, are routinely recorded, and typically, associations between genetic factors and these secondary traits are studied too. Analysing secondary phenotypes in case-control studies may lead to biased genetic effect estimates, especially when the marker tested is associated with the primary phenotype and when the primary and secondary phenotypes tested are correlated. Several methods have been proposed in the literature to overcome the problem, but they are limited to case-control studies and not directly applicable to more complex designs, such as the multiple-cases family studies. A proper secondary phenotype analysis, in this case, is complicated by the within families correlations on top of the biased sampling design. We propose a novel approach to accommodate the ascertainment process while explicitly modelling the familial relationships. Our approach pairs existing methods for mixed-effects models with the retrospective likelihood framework and uses a multivariate probit model to capture the association between the mixed type primary and secondary phenotypes. To examine the efficiency and bias of the estimates, we performed simulations under several scenarios for the association between the primary phenotype, secondary phenotype and genetic markers. We will illustrate the method by analysing the association between triglyceride levels and glucose (secondary phenotypes) and genetic markers from the Leiden Longevity Study, a multiple-cases family study that investigates longevity. © 2017 The Authors. Statistics in Medicine Published by JohnWiley & Sons Ltd.